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Abstract. We study time-dependent discrete-time quantum walks on 
the one-dimensional lattice. We compute the limit distribution of a two- 
period quantum walk defined by two orthogonal matrices. For the sym- 
metric case, the distribution is determined by one of two matrices. More- 
over, limit theorems for two special cases are presented. 

1 Introduction 

The discrete-time quantum walk (QW) was first intensively studied by 
Ambainis et al. [1]. The QW is considered as a quantum generalization 
of the classical random walk. The random walker in position x € Z — 
{ 0, ±1, ±2, . . . } at time t (6 { 0, 1, 2, . . . }) moves to x — 1 at time t + 1 
with probability p, or x + 1 with probability q(= 1 — p). In contrast, 
the evolution of the quantum walker is defined by replacing p and q with 
2x2 matrices P and Q, respectively. Note that U — P + Q is a unitary 
matrix. A main difference between the classical walk and the QW is seen 
on the particle spreading. Let a(t) be the standard deviation of the walk 
at time t. That is, a(t) = \/E(X 2 ) — E(X t ) 2 , where X t is the position of 
the quantum walker at time t and E(Y) denotes the expected value of 
Y. Then the classical case is a diffusive behavior, a(t) ~ yt, while the 
quantum case is ballistic, a(t) ~ t (see [1], for example). 
In the context of quantum computation, the QW is applied to several 
quantum algorithms. By using the quantum algorithm, we solve a prob- 
lem quadratically faster than the corresponding classical algorithm. As 
a well-known quantum search algorithm, Grover's algorithm was pre- 
sented. The algorithm solves the following problem: in a search space of 
N vertices, one can find a marked vertex. The corresponding classical 
search requires O(N) queries. However, the search needs only 0(\/~N) 
queries. As well as the Grover algorithm, the QW can also search a 
marked vertex with a quadratic speed up, see Shenvi et al. [2]. It has 
been reported that quantum walks on regular graphs (e.g., lattice, hy- 
percube, complete graph) give faster searching than classical walks. The 
Grover search algorithm can also be interpreted as a QW on complete 
graph. Decoherence is an important concept in quantum information pro- 
cessing. In fact, decoherence on QWs has been extensively investigated, 
see Kendon [3], for example. However, we should note that our results 



are not related to the decoherence in QWs. Physically, Oka et al. [4] 
pointed out that the Landau-Zener transition dynamics can be mapped 
to a QW and showed the localization of the wave functions. 
In the present paper, we consider the QW whose dynamics is determined 
by a sequence of time-dependent matrices, {Ut : t = 0, 1, . . .}. Ribeiro et 
al. [5] numerically showed that periodic sequence is ballistic, random 
sequence is diffusive, and Fibonacci sequence is sub-ballistic. Mackay 
et al. [6] and Ribeiro et al. [5] investigated some random sequences and 
reported that the probability distribution of the QW converges to a bino- 
mial distribution by averaging over many trials by numerical simulations. 
Konno [7] proved their results by using a path counting method. By com- 
paring with a position-dependent QW introduced by Wojcik et al. [8], 
Bafiuls et al. [9] discussed a dynamical localization of the corresponding 
time-dependent QW. 

In this paper, we present the weak limit theorem for the two-period 
time-dependent QW whose unitary matrix Ut is an orthogonal matrix. 
Our approach is based on the Fourier transform method introduced by 
Grimmett et al. [10]. We think that it would be difficult to calculate the 
limit distribution for the general n-period (n = 3, 4, . . .) walk. However, 
we find out a class of time-dependent QWs whose limit probability dis- 
tributions result in that of the usual (i.e., one-period) QW. As for the 
position-dependent QW, a similar result can be found in Konno [11]. 
The present paper is organized as follows. In Sect. 2, we define the time- 
dependent QW. Section 3 treats the two-period time-dependent QW. By 
using the Fourier transform, we obtain the limit distribution. Finally, in 
Sect. 4, we consider two special cases of time-dependent QWs. We show 
that the limit distribution of the walk is the same as that of the usual 
one. 



2 Time-dependent QW 

In this section we define the time-dependent QWs. Let \x) (x € 1) be in- 
finite components vector which denotes the position of the walker. Here, 
a>th component of \x) is 1 and the other is 0. Let \ip t (x)) £ C 2 be the 
amplitude of the walker in position x at time t, where C is the set of 
complex numbers. The time-dependent QW at time t is expressed by 



(i) 



To define the time evolution of the walker, we introduce a unitary matrix 



U t 



a t h 
ct d t 



(2) 



where a t , b t ,c t , d t G C and atb t c t dt =fc (t = 0, 1, . . .). Then Ut is divided 
into Pt and Q t as follows: 



a t b t 




Qt 





Ct d t 



(3) 



The evolution is determined by 

\*t+i) = £|x>®(P t |tMs + i)> + Q«hMz-i)>). (4) 

Let || \y) || 2 = (y\y). The probability that the quantum walker X t is in 
position x at time t, P(X t = x), is defined by 

P(X t = X ) = \\\Mx))\\ 2 - (5) 
Moreover, the Fourier transform |$t(fc)) £ [0, 27r)) is given by 



i* r t (fc)> = E e ^i^( ;c ))' 



with i = \/— 1. By the inverse Fourier transform, we have 

6^ ifcx 



liM*)> 



The time evolution of |<ft(fe)} is 



!#*(*)>■ 



|!Pt+i(*)> = U t (k) \* t (k)) 



(6) 



(7) 



(8) 



where U t (k) = i?(fc)t/ t and R(k) 



We should remark that 



e lk 

7?(fc) satisfies R(ki)R(k 2 ) = fl(fci + fc 2 ) and = R(-k), where * 

denotes the conjugate transposed operator. From (8), we see that 



|#t(fc)> = Ut-±{k)U t - 2 {k) ■ ■ ■ U (k) |!Pb(fc)> ■ 



(9) 



Note that, when U t — U for any £, the walk becomes a usual one-period 
walk, and \A(k)) = J7(fc)* l*b(fc))- Then the probability distribution of 
the usual walk is 



2"" Jl. 

^e ifc *tf(*)'|*b(*)> 

Z7T 



(10) 



In Sect. 4, we will use this relation. In the present paper, we take the 
initial state as 



IV>oOr)> 



(11) 



[a, p\ (x = 0) 
T [0, 0] (x^O) ' 

where |a| 2 + |/3| 2 — 1 and T is the transposed operator. We should note 
that \Mk)) = |Vo(0)>. 



3 Two-period QW 



In this section we consider the two-period QW and calculate the limit 
distribution. We assume that {Ut : t = 0, 1, . . .} is a sequence of orthog- 
onal matrices with L^s = Ho and U2S+1 = -Hi (s = 0, 1, . . .), where 



Let 



Hq 



/x(x;a) = 



ao bo 
co do 



ai bi 
ci di 



^(-|o|,|o|)(a;)i 



(12) 



(13) 



7r(l-;r 2 )\/|a| 2 -:r 2 

where Ia{x) = 1 if x € A, Ia{x) = if cc ^ A. Then we obtain the 
following main result of this paper: 



Theorem 1. 



T ' 



(14) 



where =>■ means the weak convergence (i.e., the convergence of the dis- 
tribution) and Z has the density function f(x) as follows: 

(i) 7/det(#i# ) > 0, then 



f{x) = fK{x;a e ) 
where \a^\ =min{|a |, aij} 
(ii) 7/det(i?ii? ) < 0, then 



L-<|a|»-|fl»+ ^ + a/? ) tD 



ao 



(15) 



/(«) = fi<(x;aoai) 



+ (ap + ap) bp ^ 



(16) 



If the two-period walk with det(H\Ho) > has a symmetric distribution, 
then the density of Z becomes /x(a;;a^). That is, Z is determined by 
either Ho or H\. Figure 1 (a) shows that the limit density of the two- 
period QW for ao = cos(-7r/4) and a\ — cos(-7r/6) is the same as that for 
the usual (one-period) QW for ao, since |ao| < |oi|. Similarly, Fig. 1 (b) 
shows that the limit density of the two-period QW for ao = cos(7r/4) 
and ai = cos(7r/3) is equivalent to that for the usual (one-period) QW 
for oi, since laol > ail. 
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(a) ao = cos(-7r/4), (b) ao = cos(-7r/4), 

ai = cos(7r/6) ai = cos(7r/3) 

Fig. 1. The limit density function f(x) (thick line) and the probability distribution at 
time t — 500 (thin line). 



Proof. Our approach is due to Grimmett et al. [10]. The Fourier trans- 
form becomes 

|&t(fc)) = (ffi(fc)j (i))' |*b(*)> , (17) 



(18) 



where H y (k) = R(k)H-y (7 = 0, 1). We assume 



cos U 1 sm u- 
sin 0-y — cos ( 



with 8 1 7^ zp- (n £ Z) and #0 7^ #1. For the other case, the argument 
is nearly identical to this case, so we will omit it. The two eigenvalues 
x i( k ) U = 0, 1) of Hi{k)H (k) are given by 



\j(k) = C1C2 cos2fc + S1S2 + (— 1) J iyl — (cic 2 cos2fe + S1S2) 2 , (19) 



where c 7 
Aj(fc) is 

l«j(*)> 



sin# 7 . The eigenvector |«j(fc)} corresponding to 



„ „ „2ik 

S\C2Z — C\S2 



|— C1C2 sin 2fc + (— iy \/T— (C1C2 cos 2fc + siS2) 2 | i 



(20) 



The Fourier transform !?b(fc)} is expressed by normalized eigenvectors 
\vj(k)) as follows: 



|^ (fc)) = Efe( fc )!^(fc))l^(fc)). 
3=0 

Therefore we have 

|*at(*)> = (Hi(k)H (k)y \%(k)} 



Y,^(ky{vj(k)\Mk))\v 3 (k))- 

3=0 



(21) 



(22) 



The r-th moment of X^t is 

E((X 2t ) r ) = ^> r P(X 2t = x ) 



= ^ t (fc)| i^(fc))) 

Jo 3=0 



(*)) r <«iWI*o(*)> 



(23) 



where D = i(d/dk) and (t) r = i(i - 1) x • • • x (t - r + 1). Let hj(k) = 
D\j(k)/2\j(k). Then we obtain 

E((X 2t /2tY) -> / §^^(fc)|feWl^oW)| 2 (t^oo). (24) 

Substituting ftj(fc) = x, we have 

lim £((X 2t /2i) r ) = / x r f(x)dx, 

J-\c ( \ 

where 



(25) 



f(x) = fi<(x-c t ) 



|2 2 (a/3 + a/3)si 



l- |ar-i/jr + . 



Cl 



(26) 



and |c{ | = | cos 0^1 = min{| cos#o|, | cosflil}. Since f(x) is the limit den- 
sity function, the proof is complete. □ 



4 Special cases in time-dependent QWs 

In the previous section, we have obtained the limit theorem for the two- 
period QW determined by two orthogonal matrices. For other two-period 
case and general n-period (n > 3) case, we think that it would be hard to 
get the limit theorem in a similar fashion. Here we consider two special 
cases in the time-dependent QWs and give the weak limit theorems. 



4.1 Case 1 



Let us consider the QW whose evolution is determined by the following 
unitary matrix: 



Ut 



de 



(27) 



with a, b,c,d € C. Here w t € R satisfies w t +i + Wt = m, where k\ G R 
and R is the set of real numbers. Note that ki does not depend on time. 
In this case, Wt can be written as w t = (— l) J (wo - 



+ Therefore 



the period of the QW becomes two. We should remark that 
is a unitary matrix. Then we have 



a b 
c d 



(=U) 



Theorem 2. 



Xt 
t 



Zi, 



(28) 



where Z\ has the density function fi(x) as follows: 

~o™o _i_ aab/3e~ tw ° 



h{x ) = f K{x ,a){l-\\ a \'-\f3\ 2 + ^^ 



a- 



(29) 



Proof. The essential point of this proof is that this case results in the 
usual walk. First we see that Ut can be rewritten as 



U t = 



- e iw t /2 q 




a b 




■ c ™ t /2 Q 


o e - lwt/2 




c d 




o e - lwt/2 _ 



- it 



(t) 



(30) 



From this, the Fourier transform \& t (k)) can be computed in the follow- 
ing. 

\Mk)) = {mm (^i) ur (^i) } { R{ k)R (Hta) ™ (Hta) } 

x{*(^)*(*)il(^)tf} 
x...x{fl(f)B(*)fl(f)^}fl(f)|«b(*)> 

= « (- y) W + «i/2)tf}* -R ( y) l*o(fc)> ■ 



Therefore we have 

f 2w rib /•2tt+ki/2 jr. 

ltM*)> = / ?e to |«r t (Jfe)> = / « e i(*-«i/a)« |^ t(fe _ 

JO Z7T Jk.,/2 2n 



(31) 



Kl/2)) 



R 



' 2 



'ki/2 
2tt+«;i/2 

a,i ikx / r,n \TT\t \,t,R I 



l/2 



2tt 



e lRX {R(k)Uy |!^o (fc)> , (32) 



where |<£(f(fc)} = 7? (^) |<Pb(fc - Ki/2)). Then the probability distribu- 
tion is 



P(X t = x) 



" /2 ( £T 1/2 ^ e '" {mu)t imk)) ) R (?) 



x < e 



l/2 



2tt 



2tt+ k1 /2 

V fca: [/(fc)* |<W> 



/2 



2tt 



(33) 



where U(k) = R(k)U. This implies that Case 1 can be considered as the 
usual QW with the initial state \&o(k)) = R^) \&o(k - Ki/2)) and 
the unitary matrix U. Then the initial state becomes 

\&o (fc)> = T [e wa/2 a, e- lwa/2 /3], (34) 

that is, 

\M*)) = \ Tr , , ■ (35) 

l T [0, 0] (a^O) 

Finally, by using the result in Konno [12, 13], we can obtain the desired 
limit distribution of this case. □ 



4.2 Case 2 



Next we consider the QW whose dynamics is defined by the following 
unitary matrix: 

a be iwt ' 



d 



(36) 



Here w t G R satisfies w t +i = w t + Ki, where K2 6 1 does not depend 
on t. In this case, w t can be expressed as w t = n^t + wo- Noting U t = 
R (^j-) UR (— we get a similar weak limit theorem as Case 1: 



Theorem 3. 

^ => Z 2) (37) 
where Z2 has the density function fi (x) as follows: 



f2(x) = f K (x;a)\l- I |a| 2 -|/?| 2 + 



(38) 



If Wt = 2-ntjn(n = 1,2,...), {(7t} becomes an n-period sequence. In 
particular, when n = 2 and a,b,c,d G R, {?7t} is a sequence of two- 
period orthogonal matrices. Then Theorem 3 is equivalent to Theorem 
1 (i). 



5 Conclusion and Discussion 

In the final section, we draw the conclusion and discuss our two-period 
walks. The main result of this paper (Theorem 1) implies that if det(HiH ) 
> and min{|ao|, ai|} = |oo|, then the limit distribution of the two- 
period walk is determined by Ho- On the other hand, if det(HiHo) > 
and min{|ao|, |oi|} = |oi|, or dot(HiHo) < 0, then the limit distribution 
is determined by both Ho and Hi . 

Here we discuss a physical meaning of our model. We should remark that 
the time-dependent two-period QW is equivalent to a position-dependent 



two-period QW if and only if the probability amplitude of the odd posi- 
tion in the initial state is zero. In quantum mechanics, the Kronig- Penney 
model, whose potential on a lattice is periodic, has been extensively in- 
vestigated, see Kittel [14]. A derivation from the discrete-time QW to the 
continuous-time QW, which is related to the Schrodinger equation, can 
be obtained by Strauch [15]. Therefore, one of interesting future prob- 
lems is to clarify a relation between our discrete-time two-period QW 
and the Kronig-Penney model. 
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